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Abstract. In ie O and 40 Ca an isoscalar, low-energy dipole transition (IS-LED) exhausting approximately 
4% of the isoscalar dipole (ISD) energy-weighted sum rule is experimentally known, but conspicuously 
absent from recent theoretical investigations of ISD strength. The IS-LED mode coincides with the so- 
called isospin-forbidden El transition. We report that for N — Z nuclei up to 100 Sn the fully self-consistent 
Random-Phase- Approximation with finite-range forces, phenomenological and realistic, yields a collective 
IS-LED mode, typically overestimating its excitation energy, but correctly describing its IS strength and 
electroexcitation form factor. The presence of El strength is solely due to the Coulomb interaction between 
the protons and the resulting isospin-symmetry breaking. The smallness of its value is related to the form 
of the transition density, due to translational invariance. The calculated values of El and ISD strength 
carried by the IS-LED depend on the effective interaction used. Attention is drawn to the possibility that 
in TV 7^ Z nuclei this distinct mode of IS surface vibration can develop as such or mix strongly with skin 
modes and thus influence the pygmy dipole strength as well as the ISD strength function. In general, 
theoretical models currently in use may be unfit to predict its precise position and strength, if at all its 
existence. 

PACS. 24.30.Gd Other resonances - 21.60.Jz Nuclear density functional theory and extensions - 21.30.Fe 
Forces in hadronic systems and effective interactions - 25.30.Dh Inelastic electron scattering to specific 
states 



1 Introduction 

It is an experimental observation that nuclei undergo iso- 
scalar dipole (ISD) transitions in the lhui regime of exci- 
tation energy [T]. The low-energy ISD strength function 
in N > Z nuclei has received much attention, either indi- 
rectly in the context of neutron-skin modes or in the form 
of toroidal modes [2] • Absent from recent theoretical inves- 
tigations are the IS low-energy dipole (IS-LED) modes of 
N = Z nuclei, notably 16 and 40 Ca, where the J 77 = 1~, 
T = states at 7.12 MeV and 6.95 MeV, respectively, ex- 
haust approximately 4% of the ISD energy- weighted sum 
rule (EWSR) [li] and carry little, but not negligible, El 
strength. Note, for example, that the 6.95 MeV state of 
40 Ca carries practically all pygmy dipole strength below 
10 MeV in this nucleus [5] (a weaker 1~ transition at 
5.90 MeV has been attributed to a rotational band [5]). 

The electroexcitation of the above states reveals a diffrac- 
tion minimum in the longitudinal form factors [7,8,9 . 
In the mid-seventies and for some years thereafter much 
theoretical effort was directed at accounting for their El 
strength, as a measure of isospin-symmetry violation, and 
form factors. A weak isospin mixing in the ground or 



a Email: 
darmstadt.de 



panagiota.papakonstantinou@physik.tu- 



excited state could explain the findings, though the nu- 
merical results were very sensitive to the input, including 
the isospin mixing assumed [TUlfTTlfT^lfT5] . Only a scant 
few explicit reports exist on these states within the self- 
consistent Random-Phase Approximation (RPA), mostly 
regarding their position and isospin content |14U15| . 

In this work we present self-consistent RPA calcula- 
tions of the IS-LED states in N = Z spherical, closed-shell 
nuclei. We focus on such nuclei, because of the detailed 
experimental data available for 16 and 40 Ca and be- 
cause the IS-LED states are found less fragmented. Thus 
an unambiguous comparison with the experimental spec- 
trum is possible. We use finite-range interactions, both 
phenomenological (Gogny and Brink-Bockcr) and real- 
istic or semi-realistic (unitarily transformed AV18, with 
or without a phenomenological three-body term). In self- 
consistent Hartrcc-Fock-RPA (HF-RPA), such as employed 
here, the Coulomb interaction and resulting isospin mix- 
ing is either included in both the ground and the excited 
states or ignored altogether, unlike, e.g., the valence shell 
model, where a separate treatment of the ph energies and 
wavefunctions and the ph interaction is customary. We 
will find that the IS-LED state in all studied nuclei is 
a collective lhui transition. The interactions we use tend 
to overestimate the excitation energy and some the El 
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strength of the IS-LED state, but the overall good com- 
parison with the available data suggests that the model ac- 
counts for the correct physics. The presence of El strength 
and the smallness of its value, despite the collectivity of 
the state, are elegantly explained as due to the Coulomb 
interaction and translational-invariance requirements, re- 
spectively. We attempt predictions for the IS-LED states 
in the unstable nuclei 56 Ni and 100 Sn. 

As we will discuss, there are indications that other 
RPA models do not produce any prominent IS-LED. This 
may be due to not fully self-consistent calculations. At 
any rate, no effective interaction or microscopic model, in 
general, has been tuned to describe the IS-LED properties 
and few have been tried at all. Thus the following ques- 
tion arises: can theoretical models currently in broad use 
describe correctly the low-energy ISD and El strength in 
stable, experimentally well-studied nuclei, and if not, to 
what extent can they be relied upon to describe the prop- 
erties of exotic dipole modes in, e.g., very neutron-rich 
species. Although we will not attempt an answer to the 
second part of the question here, we regard the present 
exploratory work as an initiative to start examining such 
critical matters. 

In Sec. Owe present the formalism used in this work. 
In Sec. 03 results on the 1~ response of selected nuclei are 
presented and the properties of the IS-LED are analyzed 
and compared with experimental measurements. In Sec. H] 
they are discussed in the broader context of effective in- 
teractions and of pygmy dipole strength. We summarize 
in Sec. [5j 



2 Theory 

We employ the self-consistent HF-RPA for closed-shell nu- 
clei. The HF problem is solved within a single-particle ba- 
sis spanning 13—15 harmonic-oscillator shells. The same 
effective interaction is used to construct the RPA equa- 
tions, solved within the HF basis. In particular, we employ 
a two-body Hamiltonian of the form 

H = T + y NN + V Co ui + V p , (1) 

where T is the intrinsic kinetic energy, Vnn a nucleon- 
nucleon interaction excluding the Coulomb term Vboui act- 
ing between protons and 

V p = t 3 (l + x 3 )S(r)p a (R) (2) 

is a density-dependent contact interaction (r the relative 
and R the center-of-mass position vector of the interacting 
nucleon pair). For a = 1, V p is equivalent to a three-body 
contact interaction. 

We employ various finite-range NN interactions, both 
phcnomcnological and realistic ones. We will present re- 
sults mainly with the Gogny D1S [161 parameterization 
and a unitarily-transformed AV18 realistic potential, sup- 
plemented with a phenomenological three-body contact 
term [17], which we will label here UCOM(SRG) S: 5 3 n- 
The latter is determined by transforming the S— waves 



of AV18 using the unitary correlation operator method 
(UCOM) with correlation functions determined via the 
similarity renormalization group (SRG). We also use an- 
other AV18-based potential, which we will label SRGs,53N, 
resulting from transforming the S— wave channels of AV18 
using SRG [17] . The three-body term is determined in 
both UCOM(SRG) s ,«53N and SRG S ,53N such that ground- 
state properties are well reproduced throughout the nu- 
clear chart within perturbation theory (C' 3 n = 6<3 = 
2200 and 2000 MeV fm 6 , respectively, x 3 = 1.0 in both 
cases). The flow parameters are 0.16 fm 4 and 0.10 fm 4 , 
respectively. For comparison, the pure two-body UCOM- 
transformed AV18, UCOM var , already employed in sev- 
eral studies, e.g., [T81IT5] . is adopted as well. Regarding 
the UCOM and SRG procedures and the properties of the 
various interactions we refer the reader to ref. |20| . Finally, 
for the £— closed nuclei 16 and 40 Ca we have used also 
the two-body central Brink-Boeker Bl interaction. 

The ISD response is calculated for the transition op- 
erator 

A "5 

1S D =J2( r " ~ ^ r2 ) r ^m(^i) (3) 
»=1 '* 

and the electromagnetic response using 

Z N N Z 

Oei = -j r nYi m {n n ) - -r ^2 rpYi m (ttp) (4) 

n—l p—1 

in an obvious notation, where the subscripts p and n refer 
to protons and neutrons, respectively. We calculate the 
excitation strength, B(E1 f). The above operators include 
corrections to explicitly restore translational invariance. 
However, because our calculations are fully self-consistent, 
we obtain practically the same values of strength if we use 
the uncorrected forms of these operators, 

A Z 

o I ( s ° rJ = E r i y i™(^) > ojg = 5>„w«p), (5) 

i— 1 p—1 

except of course for the spurious state, which appears at 
practically zero energy. 

Electroexcitation cross sections are calculated by using 
the proton transition density, 5p p (r), 

6p(r)=Sp p (r). (6) 

We warn against the use of a corrected isovector (IV) tran- 
sition density instead, using effective charges as in 

Z N 
Spw{r) = -^6p n {r) - —5p p {r). (7) 

The two procedures yield very different results, perhaps 
coinciding close to the photon point in self-consistent cal- 
culations. The reason is that effective charges as above 
have been derived and are only applicable for use in the 
long- wavelength limit. It is of course important to remove 
the spurious admixtures, if any, from the wavefunction or 
transition density before calculating the form factor - or 
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at least use a different prescription for testing and effec- 
tively correcting the form factor [21] . 

For the longitudinal form factors in plane-wave Born 
approximation (PWBA) we use the convention 




Note that some results reported in the literature may be 
lacking factors 1/Z or \^4it, or other. Eventually, the cross 
section divided by the Mott cross section takes the place 
of the form factor squared, in distorted- wave Born approx- 
imation (DWBA). 

3 Results 

3.1 General features of the dipole response 

In fig. [1] we show the ISD and El strength functions of the 
four N = Z nuclei 16 0, 40 Ca, 56 Ni and 100 Sn, calculated 
using the Gogny D1S and the UCOM(SRG)s,<53N interac- 
tions. The position of the giant dipole resonance (GDR) 
peak, marked by arrows in the figure, is well reproduced 
by both interactions. A feature that we find with all in- 
teractions mentioned in Sec.[3]is a strong ISD state at low 
energy, which we identify as the experimentally observed 
(in 16 and 40 Ca) IS-LED, or isospin-forbidden El transi- 
tion. It lies below the unperturbed spectrum (not shown) , 
which does not contain such strong low-lying structures, 
and it carries very little El strength. The precise position 
and strength of the IS-LED depend on the interaction. 
Strong ISD states appear throughout the spectrum. At 
high energies, they are predominantly IS and can be iden- 
tified as parts of the dipole compression mode. Strong IS 
states though appear also between the GDR and the com- 
pression mode, especially in the lighter nuclei. We also 
notice that some IS and El strength appears between the 
IS-LED and the GDR, and more so for larger A. By in- 
specting the corresponding transition densities, we find 
that some of these intermediate states are predominantly 
isovector, while others resemble rather IS proton-skin os- 
cillations. Although these states are obviously interesting 
in the context of pygmy dipole, as well as IS and toroidal, 
strength studies, we shall refrain from further analyzing 
them in the present work and focus on the IS-LED modes. 

In the following we will first present the basic proper- 
ties of the IS-LED. Then we will argue that the IS-LED is 
a collective, distinct mode of vibration. Furthermore, we 
will corroborate the correspondence we have made with 
the experimentally observed isospin-forbidden states by 
examining the electroexcitation cross section, for which 
data exist. 



3.2 Properties and collective nature of the IS-LED 

In fig.[2]we show the basic properties of the IS-LED, calcu- 
lated with the interactions Gogny D1S, UCOM(SRG)s ! 53N, 
SRGs.53N and UCOM var and how they vary with mass 



number A. For the t— closed nuclei 16 O and 40 Ca we also 
show results with the Brink-Boeker Bl potential, which 
contains no spin-orbit term. We show, in particular, the 
excitation energy, percentage of the ISD EWSR, abso- 
lute El strength, and percentage of the TRK sum car- 
ried by the IS-LED. Experimental data for 16 and 40 Ca 
are included. Regarding the toroidal nature of the IS-LED 
mode, we find that it carries approximately 15 — 30% of 
the total vortical strength corresponding to the IS con- 
vection current [23] and 4 — 12% of its energy-weighted 
sum. 

The IS-LED mode carries a non-negligible amount of 
IS strength, namely 3—13% of the energy- weighted strength 
- depending more on the interaction used and less on the 
nucleus - in agreement with existing experimental data 
for 16 and 40 Ca. The Gogny D1S interaction overesti- 
mates the excitation energy, but otherwise agrees very well 
with experiment. All AV18-based Hamiltonians give very 
similar results, indicating consistency among the differ- 
ent unitary-transformation schemes. These Hamiltonians 
overestimate the El strength. They also overestimate the 
IS EWSR and the energy, so that the absolute IS strength 
may be considered in fair agreement with experiment. 

Both Gogny and AV18-based interactions show the 
same systematics for the ISD EWSR, but not for the other 
quantities. Nevertheless, they all predict that the percent- 
age of energy-weighted El strength it carries (TRK sum) 
increases by one order of magnitude in going from 16 
to 100 Sn - which means that the absolute El strength 
increases by two orders of magnitude. 

The properties of the IS-LED do not vary smoothly 
with A. Note that the properties of another, collective 
low-energy mode, namely the 3^, as a function of A show 
a complicated pattern related to shell closures [24] , which 
would hardly show up in the subset of nuclei examined 
here. Let us also point out that in each case we have taken 
into account only the single lowest-energy and predomi- 
nantly IS eigenstate, even if a secondary IS peak is found 
nearby. The fact that 16 O and 40 Ca are I— closed magic 
nuclei, while 56 Ni and 100 Sn are not, could also play a role 
in generating the kinks in fig. [5] 

One notices the rather different behaviour of the re- 
sults obtained with the Brink-Boeker Bl interaction, in 
particular as regards the energy of the IS-LED. With van- 
ishing spin-orbit interaction and splittings the IS-LED en- 
ergy is mainly determined by the value of ITiuj, a global 
estimate of which, 41 A -1 / 3 MeV, is also shown in fig. [2] 
The above observation points to a possible role played by 
the spin-orbit coupling in determining the properties of 
the IS-LED. 

The large IS strength of the IS-LED already indicates 
a rather collective, coherent transition. We first quantify 
the collectivity of the IS-LED using a criterion proposed 
in ref. [35] • In particular, we look at the ratio N*/N p h, 
where N p h is the number of ph configurations available in 
the model space and N* the number of ph configurations 
that contribute to the state in question with an amount 
of spectroscopic strength 

S ph = \X ph \ 2 -\Y ph \ 2 >1/N ph , (9) 
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Fig. 1. ISD and El response of the Z = N nuclei 16 , 40 Ca, 56 Ni, and 100 Sn, within self-consistent HF-RPA using the 
interactions Gogny D1S and UCOM(SRG)s^3n, a transformed AV18 plus three-body term (see text). The El strength func- 
tion has been folded with a 0.5 MeV-wide Lorenzian, for visibility. Arrows mark the position of the main GDR peak (from 
photoabsorption cross section data [22]; for 56 Ni the datum for 58 Ni is given) and, for le O and 40 Ca, also the IS-LED. 



where X p h , Y p h are the RPA transition amplitudes to the 
IS-LED. N* — N p h can only be achieved for very collective 
states, where all configurations contribute with a statisti- 
cal factor 1/Nph, a really exceptional occurance. N* = 1 
might indicate an excitation generated by one dominant 
ph configuration, though not conclusively. A larger N* , 
though, should be a reasonable indicator of collectivity. 
We find, for example, for the 9.82 MeV state of 40 Ca us- 
ing the Gogny D1S interaction, N*/N ph = 11/180, to be 
compared with the value 15/180 for the first GDR peak at 
19.34 MeV. Ten out of the N* = 11 configurations are IUuj 
configurations, with a hole in the sd shell and a particle 
in the pf shell. In total there are 16 lhui configurations. 
If we solve the RPA only within the restricted Ihuj space, 
we obtain the ratios 5/16 for the IS-LED at 10.43 MeV 
and 6/16 for the strong GDR peak at 21.95 MeV. Similar 
results are obtained for the other nuclei. We conclude that 
the IS-LED is a collective valence transition. 

As pointed out and demonstrated in ref. [2B], it is im- 
portant to look not only at the collectivity of a transi- 
tion in terms of S p h contributions, but also its coher- 
ence. We find that several ph configurations, not neces- 
sarily with large S p h, contribute an appreciable magni- 



tude of (Xph — ^)/i)(p||Oisd||/i)- The majority of those 
contribute with the same sign, i.e., coherently, to the ef- 
fect of a large total ISD strength. At the same time, many 
ph configurations contribute an appreciable amount of 
(Xph — Yph)(p\\r\\h). As a rule, proton and neutron config- 
urations with otherwise the same quantum numbers con- 
tribute with the same sign. The summed positive contribu- 
tions, however, cancel out the summed negative contribu- 
tions, and thus the total amplitude of the translation op- 
erator vanishes and translational invariance is conserved. 



3.3 Structure of the IS-LED vibration 

The proton and neutron transition densities, as well as 
the isoscalar velocity fields, of the IS-LED at 9.8 MeV 
and the compression mode at 35.0 MeV in 40 Ca calculated 
with the Gogny D1S interaction are shown in fig. El We 
notice that for the IS-LED the proton and neutron tran- 
sition densities coincide much better than for the com- 
pression mode, making it an almost perfectly IS mode. 
Both transition densities are characterized by a node at 
the surface. The velocity field of the IS-LED does not 
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Fig. 2. Properties of the IS-LED mode in iB 0, 40 Ca, 56 Ni, and 100 Sn as a function of mass number. Theoretical values correspond 
to only the first eigenenergy (strongly IS) obtained with various interactions (see text), (a) Energy; experimental values from 
[27II28] . (b) Percentage of the IS EWSR; data from [SHU- ( c ) and (d) Electromagnetic excitation strength and percentage of 
TRK sum rule; data from [271128] . 



follow the compression pattern of the high-lying state, 
but involves the formation of a torus around the sur- 
face. Macroscopically speaking, the IS-LED involves the 
translation of a core (r < 2 fm) against a surface layer, 
both isoscalar and uncompressed, generating a toroidal 
surface oscillation. We find an analogous phenomenon in 
the IS OHuj quadrupole state, which is characterized by 
the same transition density as the quadrupole giant reso- 
nance, but at the same time by vortical rather than hy- 
drodynamical velocity fields (see, e.g., [29]). Let us keep 
in mind that, although collective, these low- lying dipole 
and quadrupole states carry only a fraction of the energy- 
weighted strength in their respective IS channels. As fol- 
lows from the analysis in ref . |23] , they will therefore carry 
a fair amount of vortical strength. We have also studied 
the alternative toroidal operator M tor i of ref. [30] (1st 
term on r.h.s. of eq.(25) there), essentially correspond- 
ing to out-of-phase oscillations of spin-orbit partners. We 
found that most of its strength is exhausted by a few 
states between the IS-LED and the GDR. Those same 
states carry very little ISD or B(E1) strength. The IS- 
LED mode carries a negligible amount of M tor i strength. 

For 16 and with all but the UCOM var interactions 
we find that the node of the proton transition density ap- 
pears at r = 3.1 — 3.2 fm, in excellent agreement with mea- 
sured transition densities, [3"TUT3"] . as well as a collective 
model [3lfT3] . Our results are shown in fig. [4] UCOM var 
undcrprcdicts the nuclear radius [18J ■ and therefore pro- 
duces a node at smaller r. 
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Fig. 3. Proton and neutron radial transition densities (ver- 
tical lines mark the calculated point-nucleon mean square ra- 
dius) and isoscalar velocity fields for 40 Ca and the Gogny D1S 
interaction, for the IS-LED and the compression mode. The 
distance between successive velocity vectors along each radial 
direction is 0.5 fm. 
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In general, we obtain similar results regarding the na- 
ture of the IS-LED with all interactions and for all nu- 
clei. There is a qualitative difference, however, between 
the Gogny D1S interaction, on the one hand, and the 
AV18-based potentials, on the other, leading to larger El 
strengths in the latter case. As demostrated in fig. [5] (com- 
pare UCOM(SRG)s,53N result on 40 Ca with the respective 
Gogny D1S results in fig. El, the transition densities are 
less perfectly isoscalar for UCOM(SRG)s,53N- Especially 
for larger A they show an intermediate character between 
an IS and a proton-skin oscillation. Possible reasons for the 
discrepancy between the two types of interaction include 
differences in the relative importance of the Coulomb in- 
teraction and the behaviour of the symmetry energy. In 
N = Z nuclei the Coulomb interaction and the symme- 
try energy compete against each other in the formation 
(Coulomb) or not (symmetry energy) of a proton skin. We 
find that the Gogny D1S interaction predicts for 100 Sn a 
thinner proton skin than the AV18-derived Hamiltonians, 
namely 0.82 fm or 18% of the predicted charge radius, 
vs. 20% of the predicted charge radius. Similarly, the sur- 
face dynamics and the local isospin character of the IS- 
LED will be determined by the interplay of Coulomb and 
symmetry-energy dynamics. 

By various quantitative estimates 32, 33, 34 the isospin 
T = 1 admixtures in the ground state of the T z = nu- 
clei under study increases monotonically and smoothly 
with Z. According to the present results, the forbidden 
El strength may not provide a good absolute measure 
of ground-state isospin mixing, because it does not vary 
smoothly with A. Within energy-density functional theory 
it has been observed that isospin mixing in N = Z nuclei 
does not correlate strongly with the symmetry energy or 
with the proton skin thickness [33,35 , in line with our 
elaborations. The problem of quantifying isospin admix- 
tures, however, is beyond the scope of the present work. 

From the form of the transition densities it becomes 
obvious that the small B(E1) value of the IS-LED re- 
sults from the cancellation of two large quantities: the 
r— weighted integral of r 2 Sp p (r) up to the node, i_, and 



the integral above the node, /+. As a result, deviations 
in the transition densities, such as presented in figs.[4j[5j 
can lead to disproportionally large differences in the El 
strength. By setting Vboui equal to zero we obtain a per- 
fect cancellation and vanishing El strength, but the tran- 
sition is no less collective for this reason. Its energy is 
almost the same as before, and its IS strength remains 
large. It is instructive to consider the following difference 
between the electromagnetic strength of a dipole mode 
and a mode of some other multipolarity, e.g. = 3~. 
Assume equal proton and neutron radial transition densi- 
ties throughout, 5p p (r) = 8p n {r), and therefore perfectly 
IS modes. In the octupole case r 2 Sp p (r) is surface-peaked 
and its r 3 — weighted integral, yielding its B(E3), will be 
large if its amplitude is large. If the mode is strong in 
the IS channel, it will be strong in the E3 channel. In 
the dipole case, however, the r— weighted integral, giving 
the B(E1) strength, will be zero due to the translational- 
invariance condition, 



5p p (r)r 3 dr 



[8p p {r) + Sp n (r)]r 3 dr = 0, 



even if the mode exhausts most of the IS EWSR. In both 
cases the isovector strength will be zero. The non-zero 
value of B(E1) simply means that the state is not perfectly 
isoscalar. 

In fig. IS] we show the electroexcitation longitudinal 
form factor of the IS-LED in 16 and 40 Ca within distorted- 
wave Born approximation (DWBA), to be compared with 
the experimental measurements [9,12,31.. The character- 
istic minimum is reproduced in both cases, and best by 
Gogny D1S. The form factor of the compression mode is 
also shown, for 40 Ca with the Gogny force. The differ- 
ence is obvious. Close to the photon point, and for the 
integrals weighted with a Bessel function ji(qr), rather 
than r, within PWBA, we have |7+| > |7_| for the IS- 
LED. For larger q, /_ gains relative strength and eventu- 
ally dominates, leading to a node in the form factor and 
the minimum in its absolute value. By contrast, for the 
compression mode i_ dominates already towards the long- 
wavelength limit, hence there is no minimum in the form 
factor. Thus similar, at first glance, transition densities, 
characterized by a node, lead to different, even qualita- 
tively, form factors. 

In fig. [5] we show also what happens if we calculate the 
form factor with effective charges, i.e., using the transi- 
tion density of eq. © instead of simply Sp p (r), eq. (J6]). 
The curves seriously diverge. Notice, finally, that even if 
Vboui = and Sp p (r) — 5p n (r), the form factor of the IS- 
LED state is not zero for q > 0, contrary to what would 
be obtained with the use of effective charges. 



4 Discussion 

We have found that the IS-LED mode is a collective, co- 
herent transition and that its energy may depend on the 
spin-orbit coupling. The El strength it carries is a rather 
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Fig. 5. Proton and neutron radial transition densities for the IS-LED in 100 Sn evaluated with the Gogny D1S and the 
UCOM(SRG)s i( 53n interactions and in 40 Ca with the UCOM(SRG)s,a3N- Vertical lines mark the calculated point-nucleon mean 
square radius. 
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Fig. 6. Electroexcitation form factor (absolute value) derived within DWBA as the square root of the cross section divided by 
the Mott cross section, for 16 and 40 Ca, using the interactions Gogny D1S and UCOM(SRG)s^3n- (a) Calculations for the 
IS-LED of 16 O compared with results extracted from the experimental transition density [31]. (b) Calculations for the IS-LED 
of 40 Ca compared with data [9,36 . The minimum of the reported best fit to the data, at g^ in = 0.039 fm~ 2 , is indicated by 
an arrow, (c) For the D1S interaction and 40 Ca, the result (b) for the IS-LED is compared: with the result when no Coulomb 
interaction is considered; when effective charges, eq. (0, are used; and with the form factor of the peak at 35 MeV (compression 
mode). 



delicate matter. As we saw, it is due to the Coulomb inter- 
action - for charge-symmetric Vnn _ whose role in break- 
ing the local isospin symmetry increases with Z = A/2, 
leading, for example, to the formation of a thicker proton 
skin [37] . There appears to be some correlation, within 
RPA, between the amount of IS and El strength predicted 
by different interactions. We have argued that the precise 
El strength carried by the IS-LED is determined by the 



interplay of the Coulomb interaction and the symmetry 
energy or its slope. The latter affects also the properties 
of the GDR. There is a good chance then that an RPA 
or other model reproducing correctly the properties of the 
IS-LED state and at the same time those of the GDR, 
should be able to describe satisfactorily the El strength 
and other low-energy dipole phenomena as well - a spec- 
ulation worth investigating in the future. As far as self- 
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consistent HF-RPA is concerned, the above demand is not 
trivial, as the quality of an RPA description of low-lying 
states and giant resonances depends in different ways on 
the properties of the single-particle spectrum. Effects be- 
yond RPA may be indispensible. 

A related question is whether the IS-LED mode has a 
counterpart in other nuclei, notably TV > Z ones, where 
a neutron skin may develop, that could oscillate against 
the isospin-saturated core. It cannot be ruled out that 
the core can undergo its own inner-core-vs-layer excita- 
tion, be it Pauli-suppressed due to the additional occu- 
pied neutron levels. Such a mode, if it develops, would 
likely carry little El strength, but feature prominently in 
the ISD strength function. Moreover, it could mix with 
a possible neutron-skin oscillation and influence its posi- 
tion and strength. The structure of the low-energy dipole 
spectrum would be richer than hitherto predicted. It is 
tempting to regard a neutron-skin mode as a special case 
of an IS-LED mode. To the extent that they both are co- 
herent ltiu states this is a valid classification. However, 
the proton transition density is very different in the two 
cases and an ideal neutron-skin oscillation, with a node- 
less proton transition density, would carry much more El 
strength. In any case, an electron scattering experiment 
would be able to establish the character of a given state. 
The above issues will be the subject of future work. 

Up to now we have focused on a self-consistent non- 
relativistic RPA with finite-range forces, both phenomeno- 
logical and realistic. Next we ask whether other models 
give similar results. By inspecting the literature on isospin- 
forbidden El strength we find no evidence that the con- 
ventional shell model produces 1~ collectivity in the ITwj 
regime, even though a T = dipole state with the correct 
B(E1) value can be obtained with suitable adjustments. 
Similarly, we cannot conclude on the RPA results reported 
in refs. [51155]. based on empirical single-particle states. 
We found that a simple RPA model with single-particle 
states generated by a Woods-Saxon potential and a sepa- 
rable dipole-dipole residual interaction, with adjustments 
for the energy of the spurious mode to vanish, does not 
produce Ihjui collectivity. 

Regarding other self-consistent ph models, to the best 
of our knowledge there are no reports of the IS-LED mode 
in TV = Z nuclei within modern relativistic RPA. Some 
Skyrme-RPA results are tabulated in [T5] and were first 
reported in [3"91H0] . The SkE4* Skyrme parameterization 
is used, which includes an additional three-body term. The 
position and 73(751) strengths of the T = states are fairly 
well reproduced for both 16 and 40 Ca. Energetically they 
appear almost 2 MeV below the first unperturbed ph con- 
figuration, therefore they could be collective, but no in- 
formation on their IS strenght is provided. Similarly, it is 
difficult to conclude on the Skyrme-RPA results of ref. [Tl] 
and those of an early application of relativistic field the- 
ory [4"T] . 

Systematic self-consistent RPA calculations, employ- 
ing standard Skyrme functionals, have been reported for 
the ISD strength function of Ca, Ni and Sn isotopes [32] . 
A rather strong IS-LED mode is apparently predicted for 



the N = Z isotopes, though it is not further analyzed as 
such. For 40 Ca its energy is fairly close to the observed 
IS-LED. The lowest ISD state of 56 Ni is not as prominent 
as that of 40 Ca and 100 Sn, not unlike our prediction with 
the Gogny D1S interaction, see fig.[U The isovector dipole 
EWSR below 10 MeV is overestimated in 40 Ca by a factor 
of 4.5. 

Using a fairly self-consistent RPA model (spin-dependent 
and Coulomb terms missing from the residual interaction 
but no restrictions imposed on the ph space) and the pa- 
rameterizations SkM* and BSk8 as representatives of the 
Skyrme species we calculated the ISD response of 16 O and 
40 Ca. We found two very weak states below 10 MeV which 
correspond, within 0.5 MeV ( 16 0) or 0.1 MeV ( 40 Ca), to 
unperturbed ph excitations of comparable El strength. 
Thus they can hardly be considered collective or account 
for the experimentally observed ISD strength, at least for 
40 Ca. However, by omitting the spin-orbit term of the 
residual interaction in this calculation we may have missed 
an important effect 14! . In view of the results of ref. [H] 
as well as investigations of inconsistency effects [43], the 
degree of consistency must be very important when study- 
ing the IS-LED within RPA. 

Finite-range interactions have not been extensively used 
to describe closed-shell nuclei. An early RPA application 
of the Gogny Dl interaction focused on high- lying excita- 
tions and neglected states exhausting less than 4% of the 
ISD EWSR [U] - approximately the amount of strength 
carried by the IS-LED. In ref. [2] the energy but not the 
strength of the IS-LED of 16 was calculated using Gogny 
and Skyrme forces. The importance of a correct spin-orbit 
splitting for reproducing the experimental energy of this 
state was stressed. 

It is worth mentioning that a core-layer dipole nuclear 
vibration has been studied macroscopically [45], but it is 
difficult to establish a quantitative connection with our 
work as the focus there was on the pygmy dipole strength 
of TV > Z nuclei. Last but not least, a strong ISD mode 
at the correct energy and with a similar transition density 
as we find has been tentatively reported for 40 Ca within 
the Extended Theory of Finite Fermi Systems [J5] . 

5 Summary and outlook 

In summary, we have found that the self-consistent RPA 
with finite-range interactions, such as Gogny parameter- 
izations and transformed Hamiltonians based on realistic 
potentials, predict the existence of an isoscalar, low-energy 
dipole (IS-LED) mode in spherical, closed-shell nuclei, in 
agreement with available experimental data. Here we have 
focused on TV = Z nuclei, but rather strong low-lying IS 
states have been detected also in heavier stable nuclei [T] 
3], below the alleged toroidal dipole mode. The IS-LED 
involves the oscillation of an uncompressed surface layer 
against a core. Although certain types of calculations ap- 
pear to not even predict the existence of such a mode, 
available experimental data corroborate our results. In 
particular, the strength and electroexcitation form fac- 
tor of the apparently collective low-energy ISD state of 
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16 and Ca can be simultaneously accounted for. The 
presence and small amount of El strength are related, 
respectively, to the Coulomb interaction and the form 
of the transition densities, dictated by translational in- 
variance. The excitation energy is systematically overesti- 
mated. The percentage of TRK sum carried by the IS-LED 
may grow by one order of magnitude when going from the 
lightest ( 16 0) to the heaviest nucleus studied ( 100 Sn) and 
the absolute El strength by two orders of magnitude. The 
percentage of IS strength changes more moderately. 

The precise trend of the energy, strength and fragmen- 
tation as we go to heavier N = Z depends on the interac- 
tion used. Of course, the effective Hamiltonians used are 
not tailored for the present delicate application. As pos- 
sible relevant properties we identify the symmetry energy 
and its interplay with the Coulomb potential, as well as 
the spin-orbit coupling. We have speculated that a micro- 
scopic model, not necessarily RPA, which correctly repro- 
duces the energy and IS strength of the IS-LED mode, 
along with the properties of the GDR, could lead to a re- 
liable description of the El strength and other low-energy 
dipole phenomena. Implications for the pygmy dipole strenj 
of N 7^ Z nuclei and its isospin structure will be the sub- 
ject of future work. 
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